Mode-locking in ac-driven vortex lattices with random pinning 
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We find mode-locking steps in simulated current- voltage characteristics of ac-driven vortex lattices 
with random pinning. For low frequencies there is mode-locking above a finite ac force amplitude, 
while for large frequencies there is mode-locking for any small ac force. This is correlated with the 
nature of temporal order in the different regimes in the absence of ac drive. The mode-locked state 
is a frozen solid pinned in the moving reference of frame, and the depinning from the step shows 
plastic flow and hysteresis. 



PACS numbers: 74.60. Ge, 74.40. +k, 05.70. Fh 

In 1971 A. T. Fiory observed steps in the current- 
voltage (IV) characteristics of ac-driven superconduct- 
ing thin films, analogous to the Shapiro steps found in 
Josephson junctions The steps were observed for volt- 
ages such that 27r(u)/ao = (p/q)Ct, with (v) the average 
vortex velocity, ao the triangular vortex lattice period, f2 
the frequency of the external ac-drive, and p, q integers. 
This is a particular case of mode-locking, where an inter- 
nal frequency of the system (in this case uj$ — (v)2n/ao) 
locks to a rational multiple of the external frequency. 
Several other systems with many degrees of freedom, 
such as charge density waves ||||] , spin density waves || , 
Josephson junction arrays Jq] and superconductors with 
periodic pinning arrays [ffljq], also exhibit mode-locking 
behavior experimentally. Recently, Harris et al. ]i"T[ ob- 
served the Fiory steps in YBCO and found that they van- 
ish when the vortex melting line is crossed. Historically, 
the experiment of Fiory motivated the landmark works 
of Schmid and Hauger J^] and Larkin and Ovchinikov 
p0| on the dynamics of moving vortex lattices. In Pfl0| 
it was assumed that at large velocities vortices form a 
perfect triangular lattice in which the effect of random 
pinning could be treated perturbatively. However, it is 
now clear that, instead of a perfect triangular lattice, 
there are several other dynamical phases of driven vor- 
tices: plastic flow |p^| , moving smectic, transverse mov- 
ing glass and moving Bragg glass |13H15| . These regimes 
have been observed experimentally [p.6| and in numeri- 
cal studies Jl7],[l^]. Thus, it may be of interest to study 
how the Fiory steps can arise in different moving vor- 
tex phases. Moreover, a peak in the voltage noise power 
spectrum at the "washboard" frequency loq has recently 
been observed experimentally in driven superconductors 
p9| , p0| . This is a signature of temporal order in the high- 
current steady states of moving vortices |i"5|pl| ] . We here 
present a numerical study of mode-locking for the differ- 
ent regimes of vortex velocities and its relationship with 
temporal order. 

The equation of motion of a vortex in position is: 



where r,-, = |r 7 - 



/)-X>* u i'( r <p) +p (*)' (1) 



Ti v = |r» — r p \ is the distance between the vortex i 
a pinning site at r p , r\ — $ °fe 2C 



is the distance between vortices 

and 

is the Bardeen-Stephen 

friction and F(t) = ^— [Jdc + J QC cos(Oi)] x z is the driv- 
ing force due to an alternating current J ac cos(flt) super- 
imposed to a constant current J c ic- A 2D thin film super- 
conductor of thickness d with d <C A, has an effective pen- 
etration depth A = 2A 2 jd. Since A is of the order of the 
sample size, the vortex-vortex interaction is considered 
logarithmic: U v (r) = -A v lnfr/A), with A v = <I>o/87rA 
|18| . The vortices interact with a random distribution 
of attractive pinning centers with U p (r) — —A p e~' Kr /^ , 
£ being the coherence length. Length is normalized by 
£, energy by A v , and time by r = r)£ 2 /A v , We consider 
N v vortices and N p pinning centers in a rectangular box 
of size L x x L y , and the normalized vortex density is 
n v = N v £ 2 1 L x L y = _B£ 2 /$o- Moving vortices induce a 
total electric field E = — v x z, with v = ^- v^. 

We study the response of the vortex lattice to an ex- 
ternal ac+dc force of the form F = [i^ c + F ac cos(f2£)]y 
p|| at T = 0, solving Eq. (1) for different values of F ac 
and Q. The simulations are for constant vortex density 
n v = 0.04 in a box with L x /L y = V3/2, and N v — 
64, 100, 144, 196, 256,400 (we show results for N v = 256), 
and we consider weak pinning strength of A p /A v = 0.05 
with a density of pinning centers being n p — 0.08. Pe- 
riodic boundary conditions are imposed on the simula- 
tion box and the resulting long-range interaction is de- 
termined by Ref . J2j| . The equations are integrated using 
a time step of At = O.OOIt, averages are evaluated dur- 
ing 131072 steps after 3000 steps for equilibration (when 
the total energy reaches a stationary mean value). 

Let us first review the behavior for F ac — 0. There are 
three dynamical regimes when increasing F e i c above the 
critical depinning force, F c , in the case F ac = 0: plastic 
flow for F c < Fdc < F p , smectic flow for 



1 




D 1 

G 



CM 
> 



1.000 ; 





0.90 0.92 0.94 0.96 0.98 1.00 




10 

f 



FIG. 1. Velocity- force curve around the main interference 
condition V = Q.ao/2-K for three typical drive frequencies fl. 
Each case show results for two values of amplitude F ac (the 
curves are shifted in Fdc for clarity). Insets show correspond- 
ing voltage power spectrum for F ac = and V ~ Vstep- Verti- 
cal dashed line in the spectral density indicates the washboard 
frequency, (a) fi = 0.5, F ac = 0.4 (left), F ac = 1.8 (right), 
(b) fi = 1, F ac = 0.4 (left), F ac = 3 (right), (c) fi = 2.5, 
Fac = 0.75 (left), Fac = 4 (right). 

F p < Fd c < F t , and a transverse solid for F t < Fd c (see 
lyl). The characteristic forces in our case are F c w 0.15, 
F p i=s 0.6 and F t « 1.2. Since the nature of mode- 
locking is related to the existence of temporal order at 
the washboard frequency loq we analyze the voltage noise 
Sy(f) = It Jo dt(Vy(t)-V)eKp(i2irft)\ 2 in each dynam- 
ical regime. In the insets of Fig. l(a-c) we show the 
spectral densities corresponding to each regime and we 
indicate the corresponding luq. In the inset of Fig. 1(a) 
we see that there is no temporal order at the washboard 
frequency, and only the typical broad band noise of plas- 
tic flow is observed Jl^pO]] . For the smectic flow regime, 
we see in the inset of Fig. 1(b) that there is a small and 
broad peak at ujq. Only for large forces, F > Ft, in 
the transverse solid regime well developed peaks appear 
at the washboard frequency [ p0[ and harmonics (see the 
inset of Fig. 1(c)). 

We now study the response of each of these dynamical 
regimes with velocity V = V(Fd c ) to the superimposed 
ac- force F ac cos(fit), for varying values of F ac for a given 
f2. We expect the main interference step (p = q = 1) to 
occur when V = V s tep = fla/2ir (i.e., ft = uio) if there is 
mode-locking. We therefore choose the values of f2 such 
that the expected step, V s tep = fla/2w, would correspond 
to velocities V belonging to a given dynamical regime 




FIG. 2. (a) Velocity-force curve around the main interfer- 
ence step for Q — 5 and F ac = 6. Inset shows the finite size 
dependence of the step width, (b) Time averaged quadratic 
mean displacements in the longitudinal direction {w y (t)), (O) 
points, and in the transverse direction {w x (t)), (A) points. 
The dash-dotted line indicates Oq. (c) Intensity of the Bragg 
peaks. For smectic ordering S(Gi), K y = 0: (+) points. For 
longitudinal ordering 5 , (G f 2,3), K y — 0: (*) points. 

of the limit F ac — 0. Each simulation is started at 
(v y ) ps 0.975f2a/27r with an ordered triangular lattice up 
to values such that (v y ) ss 1.025Oa/27t by slowly increas- 
ing the dc force Fd c with AFd c = 0.001. For low f2, for 
which we have plastic flow when F ac — > 0, we find that 
there are no interference steps in a wide range of F ac 
(shown in Fig. 1(a) for F ac /V ste p < 1 (left curve) and 
F ac /V s tep > 1 (right curve)). This is consistent with the 
observed lack of temporal order in the inset of Fig. 1 (a) 
and indicates that very large F ac are possibly needed in 
order to induce mode-locking. For intermediate Q, for 
which we have smectic flow when F ac — > 0, we find that 
there are no steps for small amplitudes, F ac jV ste p < 1, 
while there are steps for F ac /V s t ep > 1, as shown in 
Fig. 1(b) in the left and right curves, respectively. This 
means that the small washboard peak observed in the 
inset of Fig. 1(b) is not large enough to induce mode- 
locking steps for small F ac . However, this short-range 
temporal order can be amplified for intermediate values 
of F ac giving place to steps in this case. For high Q, cor- 
responding to a transverse solid regime when F ac — > 
we find that there are steps both for small F ac /V s t e p < 1 
and large F ac /V s t ep > 1 values of the ac amplitude, as 
we can observe in Fig. 1(c). This is in agreement with 
the F ac = spectral response 
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FIG. 3. Velocity-force curve around the depinning from 
the main step for Q = 5 and F ac = 6. The upper left in- 
sets show typical time averaged coarse-grained density of vor- 
tices seen from a system of reference moving with velocity 
Vstep = flao/2n, for a mode-locked state Fdc = 4.045 (left) 
and for an unlocked state Fdc = 4.085 (right). The lower 
right inset shows the corresponding typical quadratic mean 
longitudinal displacements for both cases. Dashed line indi- 
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observed in inset of Fig. 1(c). In this case, the temporal 
order is robust enough for mode-locking to be produced 
by small values of F ac . 

Let us now examine in detail the dynamics and the 
structural order within and in the vicinity of a mode- 
locked step, in the transverse solid case. In Fig. 2(a) 
we show a typical V — Fd c curve around the step. In 
the inset of Fig. 2(a) we show a finite size analysis of 
the step width for N v = 64, 100, 144, 196, 256, 400, where 
the error bars are due to the observed dependence of 
the width in three different realizations of disorder. We 
observe that for N v > 256 the step width tends to sat- 
urate in a size-independent value. To analyze the dy- 
namical behavior we define the quadratic mean displace- 
ments of vortices in directions parallel w y (t) and per- 
pendicular w x (t) to the external force, calculated from 
the center of mass position (X cm (t),Y cm (t)) as: w x {t) = 
-k EiMQ ~ ^(0)] 2 and w y (t) = ± - y,(0)] 2 , 

where x l {t) = x t (t) - X cm (t) y j/j(f) = y^t) - Y cm (t). In 
Fig. 2(b) we show the time average of these quantities, 
(w x (t)) and (w y (t)}, as a function of Fdc- Outside the 
step we observe that the transverse mean displacement 
(TMD) is limited (w x (t)) <C a 2 ,, while the longitudinal 
mean displacement (LMD) is unbounded (w y (t)) ^> a§. 
This corresponds to a state with only longitudinal dif- 
fusion (i.e., a transverse solid [|l8|). Noticeably, in the 
transition to the synchronization, the LMD freezes in 
a value, (w y (t)) <§; a 2 ,, while the TMD remains practi- 
cally constant. This mode-locking longitudinal freezing 
can also be observed as a dramatic decrease of the low 
frequency voltage noise in both directions. The mode- 
locked state is therefore a frozen solid. To study the 
translational order we calculate the structure factor as 
S'(k) = J2i exp[ik- ri(t)]\ 2 }. In this regime there are 
smectic order peaks of magnitude S s = S(G\) with 
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FIG. 4. Step width A fdc vs F ac 2ir/Q,a = F ac /V s tep- (a) 
Q = 1. (b) Q — 5 (A) points and £1 = 2.5 (O) points. Dashed 
line shows a fit to A\Ji (F a c/V s te P )\ and the dot-dashed line a 
fit to B|Ji(F ac /V st(!p )| 2 . 

Gi = (±27r/ao,0) and longitudinal peaks of magnitude 
Si = (S(G 2 ) + S(G 3 ))/2 with G 2 = ±27r/ao(l/2, V3/2) 
and G 3 = ±27r/a (-l/2, %/3/2). In Fig. 2(c) we plot S s 
and Si. We do not observe any important change in the 
transition to the mode-locked state. Inside the steps we 
observe jumps in Si indicating that there are different 
metastable mode-locked structures. In Fig. 3 we ana- 
lyze in detail the process of depinning from the step. We 
show a detailed view of the V — Fd c curve in the tran- 
sition from the synchronized regime. Varying Fd c back 
and forth, we observe a clear hysteresis cycle. In the 
lower right inset of Fig. 3 we show the evolution of w y {t) 
with time inside the step and outside the step. While 
w y (t) <SC eto for all t inside the step, outside there is ballis- 
tic diffusion w y (t) ~ i 2 . To visualize the spatial structure 
in the transition we define a coarse-grained vortex den- 
sity p' v (v, t) seen from a system of reference moving with 
velocity V ste p as follows: p' v (r,t) = jj- £\ 6(r - r'j(f)), 
where r'j(t) = r-j(i) — yV ste pt- We take a coarse-graining 
scale Ar < a . In the upper left insets of Fig. 3 we 
show the temporal average (p' v (r, t)) of the density, inside 
(mode-locked) and outside (mode-unlocked) the step. 
This quantity shows the stationary trajectories of vor- 
tices seen from a moving frame of reference with velocity 
V s tep- Inside the step we see that the vortices are local- 
ized, oscillating around their equilibrium positions in a 
moving lattice with velocity V s t ep - This is in agreement 
with the mode-locked frozen solid inferred from the re- 
sult shown in Fig. 2(b), which is "pinned" in the moving 
reference frame with velocity V s t ep ■ Just above the depin- 
ning from the step we see that some vortices delocalize 
following straight trajectories parallel to the force around 
"pinned" vortices, producing coexistence of mode-locked 
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and mode-unlocked channels of flow. This could be in- 
terpreted as a one-dimensional "plastic" depinning from 
the step. 

Mode-locking of the steps in the V — F dc curve can be 
characterized qualitatively by how the dc current range 
in mode-locking depends on F ac and CI, (as it was done 
for example in [^-||). In Fig. 4(a-b) we show the range 
(width) AF dc for the case F p < F dc < F t and F t < F dc re- 
spectively. The error bars and the mean values were esti- 
mated by repeating the simulation for three different dis- 
order realizations. In Fig. 4(a) we show AF dc for CI = 1 
vs F ac , which corresponds to the smectic flow regime for 
F ac 0. We see that there is mode- locking above a 
finite critical value F ac /V s tep ~ 1 (see Fig. 1(b)). For 
larger amplitudes we were not able to obtain a system- 
atic dependence on amplitude because the step widths 
depend strongly on the disorder realization in this case. 
In Fig. 4(b) we show AF dc for two frequencies SI = 2.5, 5 
vs F ac , which correspond to the transverse solid in the 
F ac = limit. We can collapse (approximately) both 
curves into a single curve if we plot AF dc vs F a c/V s te P - A 
dependence, AF dc ~ (1/C 66 )| Ji(F ac /V step )\ 2 , was found 
by Schmid and Hauger || (and also found in other elas- 
tic models) where V s t ep = clqCI /2tt and Cqq is the shear 
modulus. This result was obtained using a perturba- 
tive approach in an elastic model for the vortex lattice. 
Strikingly, our results seem to follow more closely a de- 
pendence of the form AFdc ~ A\Ji(F ac /V s t ep )\ with A 
being a constant. This dependence is the same found 
for the one-dimensional problem of an overdamped sin- 
gle Josephson junction or a particle moving in a periodic 
potential. In our case, a linear dependence of the mode- 
locking intensity with F ac would be a consequence of the 
existence of temporal order in the F ac = limit. This 
was not taken into account in the perturbative calcula- 
tion of Schmid and Hauger where mode-locking arises as 
a second order effect. 

In conclusion, it is possible to have mode-locking in 
driven vortices with random pinning for high enough fre- 
quencies, in agreement with the experiments of Fiory 
and Harris et al |ll[. The mode-locked state can be 
viewed as a frozen solid pinned in the moving frame of 
reference, and the depinning from mode-locking is plas- 
tic and hysteretic. Also, the response to an ac drive for 
different frequencies can be an interesting experimental 
probe of the dynamical regimes of driven vortices. 
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